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SUMMARY 


This investigation was performed to determine the feasibility of 
obtaining lateral frequency response and stability derivatives utilizing 
the forced oscillation technique. Tests were conducted in a North 
American Navion at one altitude, airspeed and gross weight condition. 

Both rudder and ailerons were oscillated sinusoidally with a 
simple mechanical oscillator and aircraft responses were telemetered 
and recorded on magnetic tape. A technique was developed for elimin- 
ating distortion of the input sine wave by utilizing a low frequency band 
pass filter in the reduction of flight data. Frequency responses were 
determined for D » V. and 3 for both the rudder and the 
aileron oscillations; and these were used to determine the aircraft's 
stability derivatives. Frequency responses were calculated from the 
experimentally determined derivatives and were compared with the 
actual frequency responses to illustrate the accuracy of the method 
employed to extract id stability derivatives. 

The methods employed proved entirely feasible for this aircraft 
at the single flight condition and apparently were applicable to many 
current production aircraft. The frequency responses showed little 
random error, and the responses calculated from the experimentally 


determined derivatives closely matched the actual responses. 





INTRODUCTION 


The design of modern high speed aircraft and the automatic 
control systems with which nearly all are equipped requires an 
accurate evaluation of their response characteristics. In the early 
design stages the only method available to the designer involves the 
estimation of the aircraft's stability derivatives and the solution of 
the assumed equations of motion. Because of the inherent difficulties 
involved in the analytic determination of stability derivatives, large 
errors may be present in the derived response characteristics. 
Dynamic flight testing of the prototype aircraft to accurately determine 
its stability derivatives and the coefficients of its transfer function is 
becoming increasingly important. 

Dynamic flight test data is usually obtained by one of two 
methods. The most popular of these involves the deter mination of the 
transient response of the aircraft to known control inputs. While the 
data thus obtained is, in itself, not in the best form for analysis, the 
recent development of efficient computing machines has made it possi- 
ble to convert such a transient response to frequency response form 
where it may be more readily analyzed. The second method involves 
the direct determination of the aircraft's frequency response by the 
steady state forced oscillation technique. While considerably more 


flight time is required in this method, data reduction is greatly simpli- 





fied; and good results may be obtained without the use of automatic 
computers. 

The determination of an aircraft's lateral frequency response by 
the forced oscillation technique may be accomplished by oscillating 
either of the lateral controls with a sinusoidal forcing function while 
keeping the other locked. From simultaneous recordings of the air- 
craft responses in ኃ V and / andthe control inputs on 
separate time traces, it is possible to determine amplitude ratios and 
phase angles for each response over a range of oscillation frequencies 
for the steady state condition. From these amplitude and phase spectra 
it is, intheory, possible to determine from either the aileron or rudder 
data all the aircraft's aerodynamic derivatives, except of course the 
control derivatives for the locked surface, providing the weight and 
moments of inertia of the aircraft are known. In practice, however, it 
is frequently difficult to determine all the derivatives accurately since 
some of them have little effect on the frequency response and thus are 
difficult to determine from it. 

While considerable effort has been expended on this type of testing, 
especially at the Cornell Aeronautical Laboratory, the accurate deter- 
mination of derivatives without harmonic analysis of the time traces has 
required very pure sinusoidal control inputs. This requirement has led 
to complicated instrumentation which has made the method applicable 


only to research type testing. 








This investigation was conducted to determine the applicability 
of the method to more simplified testing techniques utilizing a simple 
mechanical oscillator, telemetering equipment for data recording, and 
an ultra-low frequency band-pass filter to eliminate sine wave distor- 
tion. The flight tests were conducted in a North American Navion at 
the Forrestal Research Center, Princeton University, during the 


winter and spring of 1957. 


EQUIPMENT 


Test Airplane 


The test airplane was a North American Navion, which is an all 
metal, four place, low wing, single engine airplane. The Navion is 
equipped with a Continental E-185 engine, rated at 185 HP at 2300 
RPM and 29 inches of MP, at sea level. The wing is unswept, and the 
tail configuration conventional. The control surfaces are all of beaded 
skin construction, and have no trailing edge extensions. Trim tabs 
are located on the elevator, aileron, and rudder. 

The test airplane was modified by the removal of the two rear 
cockpit seats to provide space for the test instrumentation. The test 
airplane had a gross weight of 2870 pounds with full service, pilot and 
co-pilot, and all test equipment. 

A three view drawing of the Navion is shown in Fig. 1 and the 
general specifications of the test airplane are tabulated in Table I. 


Fig. 2 is a photograph of the test airplane. 


Power Supply 


Power available for the instrumentation in the test airplane in- 
cluded both 12 and 28 volt direct current supplied by a 28 volt engine 
driven generator and two 12 volt batteries connected in series. A 500 
VA voltage and frequency controlled inverter was installed behind the 
co-pilot's seat to provide 115 volt, 400 cycle, three phase alternating 


current power to the gyros and telemeter transmitter. 





Considerable effort was expended in obtaining an inverter which 
could be frequency controlled, since it was very important that the 
three phase power to the rate gyros be constant at the calibration 
frequency of 400 cycles. 

A circuit diagram of the aircraft 28 volt system and the inverter 


system is shown in Fig. 3. 


Aircraft Instrumentation 

Standard aircraft instruments were used to determine indicated 
airspeed and altitude. Measurement of the following quantities was re- 
quired: sideslip angle, yaw rate, and roll rate with respect to wind or 
stability axes, and aileron and rudder position angles. In choosing the 
type of components needed to measure the above quantities, an attempt 
was made to utilize instrumentation compatible with the telemetering 
equipment. For this reason, the transducers selected were of the 
potentiometer type, not only because they were capable of providing 
signals of sufficient magnitude, but also because their excitation source 
readily provided zero and full scale reference voltages, 

In view of some past history of various amounts of coupling be- 
tween the response of the aircraft and certain of the instrumentation 
components in this type of flight testing, considerable attention was 
given to choosing transducers with Satisfactory frequency response 
characteristics in order to avoid correction of response data for the 


dynamics of the recording instruments. In general, it was desirable 





that the natural frequencies of the rate gyros and the sideslip vane be 
large as compared to the expected forced oscillation frequencies of the 
aircraft. 

Sideslip angle was sensed by a light vane mounted on a boom 
approximately one maximum fuselage diameter ahead of the left wing 
tip. The moment of inertia of the vane was very low, and no gearing 
was introduced in the linkage in order to keep the natural frequency and 
damping within acceptable limits. An Electro-Mechanical Laboratory 
ultra-low torque potentiometer with .008 inch-ounces of friction torque 
was obtained for use with the sideslip vane. 

In order to reduce boom resonance, which had been experienced 
in previous tests with the Navion, the boom was stiffened by three angle 
strips equally spaced around the periphery, and running the full length 
of the boom. The sideslip vane and the boom are shown in the assem- 
bled condition in Fig. 4. 

Aileron and rudder angle were measured by 10, 000 ohm Model 
5301 single turn Helipot precision potentiometers. Due to the nature of 
the differential aileron system of the Navion, it was not considered de- 
sireable to pick off aileron angle at the aileron itself. Accordingly, the 
pickoffs were made from the control cable movement just inboard of 
the right aileron. A position close to the aileron was chosen in order 
to avoid the effects of cable stretch. The rudder deflection potentiome- 
ter was mounted on the bottom of the rudder and was actuated by a 


cable under spring tension. 





The measurement of yaw and roll rates was accomplished with 
type 3612 Giannini rate gyros. This type of gyro operates on a power 
supply of 115 volt, 400 cycle, three phase AC to the motor, and 115 volt 
single phase AC tothe heaters. The transducer is of the potentiometer 
type, with a DC resistance of 5250 ohms. These gyros were obtained 
on loan from the Instrumentation Laboratory of the U. S. Naval Air 


Test Center, NAS Patuxent River, Maryland. 


Oscillating Mechanism 


Design of the control driving mechanism was based on the follow- 
ing criteria: (a) The system must generate a sinusoidal output of con- 
trollable frequency which may be applied to either the rudder or the 
aileron control; (b) The driving link for the driven control, and the 
locking mechanism for the undriven control must each include a quick 
disconnect, not subject to binding when under load; (c) Driving and lock- 
ing links must include a trim adjustment; (d) Inter-unit rigidity and 
close tolerances are required due to the cyclical varying load; and (e) 
The unit must be compatible with the power and space limitations of the 
Navion. 

These conditions were met in the following manner. The rotation 
of a 400 cycle, 1/2 HP DC motor was converted to a pure sinusoidal 
translational motion by a scotch yoke, with suitable intermediate gear- 
ing. The driving rod to the rudder could be disconnected from the 


scotch yoke plate and bolted to the mounting plate to provide a positive 





rudder lock. The rod from the scotch yoke plate was fitted inside a 
tube from the rudder. By use of a large handled set screw the rod and 
tube could be rigidly locked at any desired trim position, and quickly 
unlocked without possibility of binding. The aileron drive was obtained 
by leading light aircraft control cable around pulleys from the horizon- 
tally translating shaft of the scotch yoke to a set of vertically translat- 
ing turnbuckles, from which a chain drove a sprocket wheel mounted 
on the control column. The sprocket was secured to the control column 
by a large handled set screw arrangement which allowed aileron trim 
adjustment and a quick disconnect. To lock the ailerons, the drive was 
opened at the turnbuckle and a cable bolted to the base plate substituted 
for the driving cable. 

The motor, gear box, and driving and locking unit was built up on 
one large section of 1/4" aluminum plate, which provided a very rigid 
base. The complete mechanism could be installed or removed as a unit, 
the plate being bolted to extrusions riveted to the siffeners under the 
deck on the starboard side of the cockpit. In-flight shifting from the 
aileron drive to the rudder drive involved about five minutes. 

Control of the driving motor was obtained by the use of an ampli- 
dyne mounted behind the co-pilot. Through control of the field of the 
amplidyne, whose output provided the armature current to the driving 
motor, any speed from zero to maximum rated, including reversing, 


could be obtained. Indication of the operating frequency was not directly 





available in the cockpit, but could be obtained approximately from a 
calibrated voltmeter across the amplidyne signal circuit. At the lower 
frequencies, the torque (proportional to the armature current) became 
marginal, and was evidenced as a distortion to the sine wave of the 
driven control as the aerodynamic load varied with control displacement. 
A worm and gear was incorporated in the gear box to obtain irreversi- 
bility between the motor and the control surfaces, but the variation of 
the frictional load on the worm was sufficient to cause this low frequency 
distortion. The aircraft instrumentation is shown in Fig. 6 anda close 


up of the oscillating mechanism is shown in Fig. 7. 


Signal Circuit De sign 


The telemeter transmitter was designed to accept voltages from 
the transducers from zero to five volts. It was thus desirable for the 
full scale voltages on the potentiometers to be maintained at five volts. 
Variations in this full scale value were permissable, providing that the 
variation on all the potentiometers was the same, since the telemeter 
corrected all voltages to a percentage of the full scale reference voltage. 
The potentiometer used in the sideslip vane was a 360 degree single 
turn model and since only 40 degrees of variation in sideslip angle were 
expected it was necessary to apply 45 volts excitation to this potentiome- 
ter in order that its output might have a five volt variation. The use of 
a 45 volt DC signal circuit supply was thus dictated. It was also desira- 


ble to maintain the signal voltage to the other potentiometers at five volts 
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and introduce a biasing circuit to provide a variation of sideslip output 
voltage from zero to five volts rather than a five volt variation about 
some intermediate voltage. 

The signal circuit design based on the preceding requirements is 
shown in Fig. 5. The essential element of the circuit is the bias sup- 
ply voltage potentiometer, pot "A", With the sideslip vane against the 
Stop in one direction, pot "A" was adjusted so that the output voltage 
of the sideslip pot was zero. The vane was then deflected to the oppo- 
site stop, and the output voltage of the sideslip pot was recorded. 

This value was slightly less than five volts. Pot "B" was then adjusted 
to bring the full scale reference voltage, and thus also the voltage in 
the remaining signal circuits in agreement with the full scale sideslip 
voltage. The procedure was repeated two or three times as necessary, 
in order to nullify any loading on the system. Once the final adjustment 
was made, the zero to full scale voltage for each quantity to be meas- 
ured was exactly the same, and small variations in signal voltage pro- 
duced the same variation for all circuits. Pots "A" and "B" were 
locked and required no further adjustment for the remainder of the 
project. 

During the initial flight and transmission tests, considerable AC 
ripple was found to enter the DC signal voltages, very likely from the 
400 cycle inverter. Condensers "C" and "D" in Fig. 5 were installed, 


and were found to reduce the ripple to a negligible value. 





Telemeter System 


The telemeter unit was of the pulse width type and was designed 
and built by the Applied Science Corporation of Princeton (ASCOP). 
This unit was the key element in the instrumentation of this project. 
The basic functional components of this pulse width system are trans- 
ducers (all of the potentiometer type); a time division multiplexer, or 
eommutator; an electronic coder, or keyer to put the intelligence into 
pulse form; an RF transmitter; and ground station equipment for re- 
ceiving, decoding, and recording the data in usable form. 

The airborne component of the system is designated the ASCOP 
D Series PW Multicoder, and the ground station equipment the ASCOP 
M Series PW Ground Station. 

A unique feature and very great advantage of pulse width data 
systems is that, due to the multiplexing or time sharing operation, all 
data channels traverse the same path during transmission. Since zero 
and full scale references are transmitted with the data, this results in 
a calibration signal, effective for all channels, being transmitted 20 
times each second. These calibration signals allow the removal of 
almost all of the system error, with the exception of errors introduced 
in the instrument pickups and the mechanical switch, and they make 
possible the use of this telemetering system in tests of long duration 


and varying environmental conditions. 





1ረ 


Magnetic Tape Recorder 
An Ampex Model 309C dual track three speed tape recorder was 


an integral component of the ground station equipment, and provided a 
convenient means for recording, storage, and playback of the standard 
pulse width modulated data. Provisions were also included for placing 
voice transmissions directly on the magnetic tape from either or both 


the aircraft and the ground station. 


Data Recorder 

Qualitative analysis of the flight data during the actual flights and 
quantitative analysis of the recorded signals was made on a Sanborn 
Model 154-100D four channel recorder. This instrument had a very 
high natural frequency (42 cps), and the low frequency response was 
flat to zero cps. Different paper speeds and variable gain were possible 
and permitted expansion of the time trace and amplification of the signal 


to nearly any desired value. 


In order to study the response of the aircraft to a true sinusoidal 
driving function, all forcing functions and the corresponding responses 
were filtered to remove all but the fundamental frequency. The filter 
employed was the model 330-A KROHN-HITE Ultra-Low Frequency 
Band-Pass Filter. This is an adjustable filter whose gain is unity in 


the pass band and drops outside the pass band at the rate of 24 db per 
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octave. The filter combines a high pass filter and a low pass filter, 
each of which may be continuously tuned from 0.02 to 2,000 cps. 
Since the filter will handle an input signal of only up to 10 volts rms 
without overloading, a voltage divider was applied to the telemeter out- 
put, which had a full scale range of up to 100 volts. For this applica- 
tion the filter was used as a narrow pass band with the high and low 
dials set at the same frequency. Used inthis manner, the filter did 
not alter the relative amplitude or phase relationships between the 
input and the output at a specific frequency. 


A photograph of the ground instrumentation is shownin Fig. 8. 


Instrument Calibration 

Six items required calibration in the instrumentation system. 
These were the two rate gyros, the two control deflection potentiome- 
ters, the sideslip vane potentiometer, and the signal voltage to the 
amplidyne which was roughly a measure of the frequency of control 
oscillation. The aircraft altitude and airspeed indicating systems were 
not calibrated because previous calibrations indicated that the errors 
here were negligible for the conditions of this investigation. 

The two rate gyros were calibrated on a pendulum which consisted 
of a 3/4" diameter aluminum shaft and a 12" diameter instrument pan. 
The pendulum, the overall length of which was 121.75" was pinned to 
the overhead truss structure in the Forrestal hangar. A scale was 


mounted with brackets on a small table in such a way that displacement 
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of the pan from its rest position could be measured. The gyros were 
mounted one at a time on a wooden base plate which was clamped to the 
pendulum pan, and their potentiometer outputs were connectedto a 
Brush Recorder, The pendulum was displaced from its rest position 
and allowed to swing freely. Since there was an excessive amount of 
friction in the system, swing amplitudes were measured on two succes- 
sive swings in the same direction and averaged. The two maximum 
voltages recorded during the swing were also averaged and the maximum 
rate attained during the swing was determined from the following 


relation: 


where P was the natural period of the pendulum whose measured value 
from the trace on the Brush recorder was 3.33 sec., 1 was the length 
of the pendulum, and d was the average maximum displacement of the 
pan from its rest position. The pendulum was allowed to swing until 
once again at rest with readings taken, as described above, over the 
entire range of gyro rates. The angular rate gradient of both gyros 
was found to be linear and equal to .610 degrees/sec. per percent full 
scale potentiometer voltage. The gyro calibration curves are shown in 
Fig. 9. It was found that the gyros were very sensitive to the frequency 
of the 400 cycle supply voltage and as a result great care was taken 
during calibration and also during the actual flight test program to 


maintain the frequency of the supply at exactly 400 cycles per second. 
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The control deflection potentiometers were calibrated by measur- 
ing the resistance between the output terminals over the full range of 
control deflections. The actual control deflections were measured with 
a vernier protractor capable of accuracy to 0.1 degree. The rudder de- 
flection gradient was found to be linear and equal to .405 degrees per 
percent full scale potentiometer voltage. The Navion aircraft is 
equipped with differential ailerons and it was expected that the aileron 
calibration would therefore be non-linear. However, using average 
aileron angle, it was found that over the small deflection range of +10 
to -10 degrees desired during the flight test program the aileron deflec- 
tion gradient was linear and equal to .1933 degrees per percent full 


scale potentiometer voltage. 


The sideslip vane potentiometer was calibrated in the same man- 
ner as the control deflection potentiometers, using a protractor and 
measuring the resistance between the output terminals with an ohm- 
meter. The calibration curve is shown in Fig. 10. The sideslip 
gradient was found to be linear and equal to .468 degrees per percent 


full scale potentiometer voltage. 


The telemeter equipment corrected its output to percent full scale 
voltage and since full scale voltage on ali the potentiometers varied by 
the same increment as the reference voltage for small fluctuations of 
the latter value, it was unnecessary to calibrate the potentiometers 


through the telemeter equipment. 


The speed control on the oscillating device was calibrated roughly 


in flight by timing the frequency of control oscillations and recording 
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the voltmeter reading across the control field of the amplidyne genera- 
tor. This speed control was used only to obtain order of magnitude 
settings, the actual frequency of any particular oscillation being deter - 
mined from the trace on the Sanborn Recorder. 

Careful consideration was given to the possible errors which 
might be introduced in the data by the dynamics of the recording instru- 
ments. The natural frequency of the rate gyros was 17 cycles per 
second and their damping ratio was between .45 and .7. The natural 
frequency and damping ratio for the sideslip vane were experimentally 
determined in flight to be 15.44 and 0.081 respectively. Since the 
natural frequency of the aircraft was approximately 0.4 cycles per 
second and the maximum forced oscillation frequency was not expected 
to exceed one cycle per second, the attenuation of the signals due to the 
dynamics of the instruments was effectively zero. While a small error 
in phase angle measurement was introduced at the higher frequencies, 
it was considerably smaller than the reading accuracy with which phase 
angle could be determined from the data, and thus was not considered 


in the reduction of flight test data. 
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FLIGHT TEST PROCEDURE 


In this investigation, frequency response was measured at the 
single flight condition of 5000 ft. density altitude, 129.6 mph true air- 
speed, and 2870 lbs. gross weight. Each test was commenced at 
steady level flight with the amplidyne control set at zero speed, the 
oscillating mechanism exactly centered to its mid position, and the 
lateral controls locked in this trim position. The run was then initi- 
ated by setting the speed control and allowing steady state conditions 
to become established. The necessity for driving the lateral controls 
about the exact trim position was so great that the lock-on system to 
a moving oscillator, as employed in a somewhat similar longitudinal 
study (Ref. 1) was not feasible. In spite of the great care exercised 
to center the oscillations about a straight and level flight path, prac- 
tically every run was ultimately terminated by the excitation of the 
spiral mode. At this time the rapid unlocking features in the mechan- 
ical system design were most appreciated. 

In setting the speed control, an approximation to the desired 
frequency could be obtained by the use of a voltmeter across the ampli- 
dyne control field calibrated to read frequency. The relatively short 
duration of each run did not permit fine adjustments to obtain a specific 
frequency, so several runs around a desired frequency (primarily near 


the dutch roll) were employed. In addition, it was found that the low 
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frequency runs were best started by going rapidly to a high frequency, 
and then coming down to the desired frequency. Otherwise, the 
response to the slow initial deflection drove the aircraft so far into the 
spiral divergence that the subsequent control reversal was insufficient 
to establish an oscillation about a straight flight path. Although it had 
been expected that all controls except the driven one would be locked 
during a run, it was found that slight elevator movements were required 
to maintain a constant airspeed throughout the run, particularly during 


oscillations of the rudder. 
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RESULTS 


Data Reduction 

The quantities measured during the flight test program were da 
عو‎ 4. , 5 , p , and 7 . These four quantities were teleme- 
tered to the ground station, recorded on magnetic tape, and simulta- 
neously plotted on the Sanborn Recorder, This raw data had noise and 
distortion superimposed on the fundamental sinusoidal driving frequency. 
An example of such a record is shown in Fig. ll. In order to study the 
response of the aircraft to only the fundamental of the driving frequency, 
a KROHN-HITE Ultra-Low Frequency Band-Pass Filter was used to 
filter both the driving function and the aircraft responses. From the 
original record of the raw data the fundamental frequency for a particu- 
lar run could be determined with sufficient accuracy to establish the 
high and low settings for the band pass filter. With the filter thus set 
to pass only the specific frequency, all four quantities were successively 
played from the magnetic tape through the telemeter receiver (convert- 
ing the pulse width coding to a varying voltage output), thence through 
the filter to the Sanborn recorder. As each quantity was being filtered 
and recorded, an unfiltered trace of the forcing function was also 
placed on this new record to act as a time reference. The results of 
this time and recorder paper consuming process were four sinusoidal 
records of the measured quantities, each with a common unfiltered time 


reference. Fig. 12 shows an example of such a filtered record with the 
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unfiltered forcing functions matched for each response. From these 
matched responses the amplitude ratios and phase angles relative to 
the driving quantity could be determined. Since the settings on the 
filter and the recorder were not altered during the filtering of a given 
frequency, any amplitude attentuation or phase shifting was applied 
equally to all four quantities. The gain of the recorder was initially 
set for each frequency to accommodate the quantity with the greatest 
amplitude, such that this quantity would utilize almost full scale on a 
recorder paper channel. By using this technique, and increasing the 
paper speed for the higher frequencies, a rather accurate determina- 
tion of the amplitude ratios and phase angles was possible. 

The quantities @ and WY were determined from the records 
of ይ and D by integration. This was readily performed by divid- 
ing by w and subtracting 90° from the phase angle, since the re- 
sponse was sinusoidal. The magnitude of each response was divided 
by the magnitude of the forcing function and the appropriate scale fac- 
tors were applied to give the desired amplitude ratios. Since the sideslip 
vane was not located at the c.g. of the aircraft an additional correction 
tothe § response was required. This correction was a function of 
the distance of the sideslip vane from the c.g. of the aircraft and is 
given by the following relation: AG = — .027¢ Y 
The final frequency responses were tabulated in Tables II and III and 


were plotted in Figs. 13-18 and labeled "experimental." 
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Also plotted on the appropriate graphs in Figs. 13-18 were fre- 
quency response curves denoted as "theoretical." These curves were 
the result of an initial study of the Navion airplane based on lateral 
stability derivatives determined from previous tests, such as described 
in Ref. 2, or from theoretically determined derivatives. These initial 
"best guess" derivatives are tabulated in Table IV, and were used as a 
rough check on the results, as well as the basis for the initial system 


design. 


Determination of Derivatives 


Ref. 3 develops a method for extracting stability derivatives from 
frequency response data for the lateral case. The method involves the 
application of the vector technique to the three linear equations of lateral 
motion to separate them into six equations. The vector responses are 
resolved parallel and perpendicular to the control deflection vectors. 
The perpendicular resolution eliminates the control deflection deriva- 
tives in one of the two equations derived from each basic equation and 
thus permits solution of the aerodynamic derivatives separately. The 
development of these equations is shown in Appendix A. The method of 
least squares is then applied to the data for each of these equations to 
yield lateral stability derivatives. 

Although the method of least squares is probably the best way to 


eliminate random error, it is not readily applicable to hand calculation 
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techniques. Automatic computing facilities were not available here, 
and since time was limited, it was decided to use the simpler method 
of averages in conjunction with the vector equations developed in Ref. 3 
to determine the lateral stability derivatives of the Navion. The K 
notation of Ref. 3 was modified to a starred derivative notation for 
greater clarity. This notation is defined in Appendix A. The method 
of averages, as described in detail in Ref. 4, is also outlined in 
Appendix A. 

The components of the responses along and perpendicular to the 
control deflection axis are defined in the form A+B with the appropri- 
ate subscripts. Formulas for the determination of the responses and 
their derivatives in this form are given in Appendix A. 

The two equations developed from the basic side force equation 
by the vector technique are equations (13) and (14) in Appendix A. 
The equation describing the relation between the components perpen- 


dicular to the control deflection axis was: 


* 
03) Cyg - C" 3$ = -w (Ag+Ay) 


Since the only unknown was owe it was possible, in theory, to deter- 
mine its value from the responses to either an aileron or rudder forcing 
function at a single frequency. Since the responses in yaw and sideslip 
were greater for a rudder input, it was felt that the rudder data would 


give a better solution for Cyg . In order to eliminate random error 
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the responses at 12 selected frequencies were taken from the faired 
rudder curves. The summation of the 12 equations determined from 
these 12 points gave one average equation, from which the value of 
2 was determined. This gave a value of Cyg of -0.591. The same 
technique using aileron data resulted in a positive value for Cyg ; 
which was, of course, impossible. Ref. 3 mentions that the accurate 
determination of Cy g may be difficult from this equation, since it 
involves the sum of the yaw and sideslip components, which are nearly 
equal but opposite in sign. The authors of Ref. 3 recommend obtaining 
lateral acceleration data and using a modified side force equation to 
avoid this condition. The large responses in yaw and sideslip to rudder 
inputs tended to minimize this condition, but the determined value of 
Cy 4 was nevertheless subject to considerable error. 

The equation describing the relation between the components 
parallel to the control deflection axis and therefore including the deriva- 


۱ * ۱ 
tive, Cys, , was: 


* % 
aa) Cyg Ag- C? Ap = Cys, + w(Ba+ By) 


Using the previously determined value for C, $ in an average equation 
for the 12 selected points yielded a value for Cy > This gavea 

value for Cys. of 0.206. The accuracy of Cy¿_ was clearly depend- 
ent on the correctness of Cy 4 , and was, in addition, subject to small 
differences between yaw and sideslip components. Cy f, was therefore 


of doubtful reliability. 
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The vector technique gave equations (15) and (16) of Appendix A 
from the basic roll equation. The relation between the components 


perpendicular to the control deflection axis was: 
> 7ታ 5 
08) Ceg Ba t Cu Ag e 15$ ^8 - Cap wA = بت‎ BS 


Since L., was very small compared to Ix the term E was assumed 
to be zero. Solution for the three unknowns required three equations. 
The 12 frequency points from the aileron data were now grouped to in- 
clude four frequencies below the dutch roll, four near the dutch roll, 
and four above the dutch roll, and average equations were determined 
for each group by adding the appropriate coefficients. Solution of these 
equations for Cle” C1,” and ent were converted to Cle ; 

E. and C1, which were -0.0505, -0.436, and 0.0390, respectively. 
In the solution of these equations it was found that an accurate determin- 
ation of Cl, was difficult, since the solution involved a small differ - 
ence between large quantities; and small errors were magnified many 
times. An attempt was made to carry greater accuracy in the solution 
than the data actually warranted, but of course this procedure, while 
required to give any solution at all, could not give an accurate solution 
for Cj,. This condition was not surprising, since it is well known 
that the cross derivatives are difficult to determine from this sort of 
data. 


Solution of the same equation using a grouping of points from the 
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rudder data should have given the same derivatives as previously 
determined from the aileron data. When this solution was attempted, 
however, it was found that the simultaneous equations were ill-condi- 
tioned. This arose from the fact that £ and & varied in nearly 
the same manner over the entire frequency range, both in phase angle 
and magnitude. Additional frequency points were taken and various 
groupings were tried in an attempt to eliminate this difficulty, but all 
were unsuccessful. Ref. 3 suggests a possible method of avoiding the 
difficulty by eliminating the cross derivative coefficient by resolution 
of the vectors along the 0 axis, instead of along and perpendicular 
to the control deflection axis. This was attempted, but solution of the 
modified equation was still impossible. 

It was felt that the rudder data was slightly better than the aileron 
data, due to refinements in the flight test technique as the program pro- 
gressed, and it was therefore highly desirable to include the rudder 
data in the solution for the roll derivatives. The solution of the aileron 
equations gave a value of C1,” which was felt to be quite accurate, 
since it was in close agreement with theory and previous experimental 
results. This value was assumed known in the rudder equations for the 
rolling motion, and the two remaining derivatives Gy," and Gua 
were determined by taking two groups of six frequency points each and 
arriving at two equations in the two desired unknowns. The values de- 


termined for “14 and Cj, were -0.0515 and 0.119 respectively. 
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It was interesting to note that the value of Clg determined in this 
manner was nearly the same as that determined from the aileron data, 
and the doubtful derivative, Cl], was nearer the theoretically deter - 
mined value. The derivatives determined in this manner were arbi- 
trarily chosen as the final values, since the authors felt that they repre- 
sented the best possible combination of all the experimental data. 

The relation between the components parallel to the control 


deflection axis was: 


+ 3 » 
(16) Cag Ag - Cah uw Bq + u eo Apt Cap v Bx, = Ce, = wrAd 


Using the previously determined values of C],*, Cj,*, and C1," 
and utilizing one average equation based on all 12 frequency points 
from the rudder data permitted solution for Cj, *, which gave a 
Cy [o of 0.0128. Similarly, solution of the same equation using the 
aileron data gave Gites *, which resulted in a value for C14, of 
-0.108. 

Expressing the basic yaw equation in vector form resulted in 
equations (17) and (18) of Appendix A. For the components perpen- 


dicular to the control deflection axis the relation was: 


* * 1 * 
E 


Because oí the basic similarity between the equations for the yawing 
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and rolling motions, the unknowns in equations (15) and (17) have the 
same numerical coefficients, with only the constant term differing. 
Because of this, rapid solution of the three average equations for the 
yawing motion was possible using the same grouping of frequency 

points from the responses to aileron oscillations. This gave values of 
Che , “እይ , and Car of 0.0771, -0.129, and -0.0841, respectively. 
The same difficulty existed here in the determination oí Cp, as had 
previously existed in the determination oí C],- 

Because the numerical values of the coefficients of equation (17) 
were the same as those of equation (15), the three simultaneous equa- 
tions for the yawing motion derived from the rudder data were also ill- 
conditioned for solution. The method applied to solve equation (15) 
using an apparently accurate value of C1, could not be applied in this 
caSe, because both Cua and Cup as determined from the aileron 
data were felt to be doubtful, since they were not in agreement with 
either theory or the results contained in Ref. 5. Inthis case final 
solution was effected by solution of a set of simultaneous equations, 
two of which were determined from the aileron data at frequencies near 
and above the dutch roll, and the third representing the average of all 
the rudder data. This approach appeared legitimate since each equa- 
tion in the set was independent of the type of forcing function. This 
particular choice of equations resulted from the opinion of the authors 


that the aileron data was better at higher frequencies and that inclusion 
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of the generally more reliable rudder data was desirable. It should be 
emphasized that this particular choice was peculiar to this experiment. 
Solution of these equations gave values for Ong : Cp , and Cr, 
of 0.0895, -0.106, and -0.0908 respectively. 

The equation relating the components parallel to the control 
deflection axis for the yawing motion was: 

(18) -G. As + وقد وی‎ + EE utag - Cu wBy - ርእ... = መ”ዶሌጭ 
Using the previously determined values of 2 ۱ ek , and c 
and utilizing all 12 frequency points from the rudder data permitted 
solution for C, ( *, which gave a Cay, of -0.0673. Similarly solu- 
tion of the same equation using the aileron data gave Cmn, ማው which 
resulted ina value for Cy , of 0. 00376. 

The final derivatives determined in this investigation are tabulated 
in Table IV. This table also contains the original theoretically deter - 
mined derivatives, and the derivatives determined for the same aircraft 
by an analog computer solution as outlined in Ref. 5. 

In order to check the accuracy of the method used to determine the 
aircraft's stability derivatives, the transfer function coefficients based 
on them were determined as outlined in Appendix A, and from these the 
frequency responses were calculated. The experimental transfer func - 


tions are tabulated in Table V, and the frequency responses based on 


them are shown in Fig. 19-24, 
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DISCUSSION 


Dynamic flight testing is becoming increasingly important with the 
advent of high speed high altitude aircraft because the aircraft designer 
is more concerned than ever before with the accurate determination of 
stability derivatives and the associated transfer functions. Although 
transient methods are the most popular at the present time, it was felt 
that if the instrumentation requirements of steady state dynamic flight 
testing could be reduced and the necessity for absolutely sinusoidal con- 
trol inputs without the added requirement of harmonic analysis could be 
eliminated the forced oscillation technique for determining frequency 
response and thus stability derivatives and transfer function coefficients 
could still be used to advantage under certain circumstances. It was the 
purpose of this investigation to examine the possibility of obtaining ade- 
quate frequency response with relatively simple instrumentation, and if 
possible to eliminate the necessity of harmonic analysis due to distorted 
sinusoidal inputs. 

Since investigation of the methods involved was of primary con- 
cern here, flight tests were conducted at only one altitude, airspeed, 
and gross weight condition, although it should be emphasized that in any 
actual flight test program to determine derivatives, all of these condi- 
tions should be varied over the operating range of the aircraft. 

Before the design of the instrumentation system could be accom- 


plished a prediction of the results based on theoretically determined 
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derivatives and control gradients was required. Obviously the experi- 
mental data had to include frequencies on either side of the dutch roll 
with frequencies at the dutch roll being of primary interest. Previous 
investigations of the lateral dynamic stability of the Navion indicated 
that the period of the dutch roll was between 2 and 3 seconds. An 
analysis based on the best available derivatives, some calculated theo- 
retically from Ref. 6 and others taken from Ref. 2 and 7 was performed. 
This indicated that the period of the dutch roll was about 2.4 seconds. 
Theoretical frequency responses to both rudder and aileron forcing 
functions were calculated and these served to design the instrumentation 
system. 

It was decided to design the oscillator to cover frequencies be- 
tween 0.5 and 10 radians per second. The upper limit proved higher 
than necessary and better results would have been obtained if the design 
had limited the upper frequency to 5 radians per second. This upper 
limit was found to be well above the dutch roll frequency and oscilla- 
tions at higher frequencies produced only negligible responses. Exam- 
ination of the control deflections necessary to produce recordable 
responses without exceeding the linear limit of the equations and the 
control gradients led to the decision to use a. à H. P. motor and an 
amplidyne to control its armature current. Early experimentation with 
a mechanical variable reduction gear indicated that this method was 


unsatisfactor y. 
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Ref. 1 recommends a careful preliminary study of the dynamics 
of the recording instruments since, if their natural frequency is not 
much larger than the natural frequency of the aircraft, large errors 
arise which must be corrected for. This was done, and it was found 
that the natural frequency of the rate gyros and sideslip vane were such 
that the attenuation of their signal outputs was effectively zero and that 
the phase shift introduced was less than 0.5%, Since other errors in 
phase angle measurement were expected to be of considerably greater 
magnitude, no attempt was made to correct for the dynamic response 
of the recording instruments. 

The instrumentation system based on the results of the above cal- 
culations proved completely adequate with the exception of the oscillator. 
It was found that since the output torque was directly proportional to the 
armature current of the motor the torque at low frequencies and thus 
low armature currents was insufficient to produce anything approaching 
a pure sinusoidal control motion. Consideration was given to changing 
the reduction gear ratio and thus providing more torque at low frequen- 
cies, but this was discarded because of time requirements. The gyros 
were properly utilized, and their limitations were never exceeded; and 
the original estimates of the control deflections required to produce the 
desired responses proved accurate, making major modification of the 
original instrumentation unnecessary. 


Several difficulties were encountered during the flight test program. 
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The initial flights were conducted with the rudder locked and with the 
ailerons driven by the oscillator. At oscillation frequencies well 
above the dutch roll the steady state condition was rapidly reached and 
satisfactory data could be obtained before the spiral mode was excited. 
At frequencies near the dutch roll the spiral mode was frequently ex- 
cited before the steady state condition was reached necessitating ter- 
mination of the run before any data could be obtained. At frequencies 
below the dutch roll the initial control deflection was sufficient to ex- 
cite the spiral mode and no data at all could be obtained. After some 
experimentation it was found that most of this difficulty could be elim- 
inated by starting the oscillations at a high frequency and then reducing 
the frequency to the desired value, and this procedure was adopted in 
all subsequent flight tests both for aileron and rudder oscillations. 
During the runs involving aileron oscillations there was very little indi- 
cation to the pilots when the dutch roll frequency was being approached, 
and it was not until the data was reduced that it could be seen that the 
responses peaked at this frequency. During the runs involving rudder 
oscillations, however, the responses became quite violent at the dutch 
roll and apparently very sensitive to airspeed. During the first rudder 
run at the dutch roll frequency the airspeed control was poor, and as 
speed decreased by as little as 6 mph the responses decreased in a 
very noticeable manner. As speed was increased to the desired value 


by retrimming the elevator, the resonant condition was again reached 
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and the responses became large. A similar variation was noted with an 
increase in airspeed over the desired value, but the effect was not as 
pronounced. The difference between the effects was explained by refer- 
ence to the final response curves which show a rather flat top at the 
dutch roll frequency. A decrease in velocity corresponds to an increase 
in 7 and this in turn corresponds to a decrease in the dutch roll fre- 
quency. Since the actual frequency of the oscillations was maintained 
constant, the operating point at the decreased velocity was actually 
above the dutch roll. Since the experimentally determined response 
curves fall off more rapidly above the dutch roll than below it, the con- 
dition was satisfactorily explained. Because of the violence of the 
responses near the dutch roll and their sensitivity to airspeed great 
care was necessary to keep the speed at its desired value and runs 
where speed control was poor were discarded. 

Ref. 1 discusses some difficulty in the design of satisfactory con- 
trol locking and trim devices for the oscillation mechanism. The orig- 
inal design accomplished here utilizing the recommendations of Ref. 1 
proved highly successful and no difficulty was encountered either in 
obtaining proper trim or in releasing the oscillator or the locked control 
surface when the attitude of the aircraft made this necessary. 

While it was felt in the original design that a measure of the con- 
trol field voltage to the amplidyne would serve as a satisfactory indica- 


tion of the frequency of oscillation, it was found that this assumption 
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was not wholly justified. A rough determination of frequency could be 
made from this voltage reading, but because of residual magnetism in 
the control field of the amplidyne, the desired accuracy could not be 
attained. Precise frequency determination for data analysis was made 
from the recorded traces but since a fairly accurate knowledge of fre- 
quency in the cockpit during the flight testing was desirable it is recom- 
mended that in subsequent tests of this nature a tachometer generator 
be incorporated on the driving motor. 

In the analysis of the data at low frequencies it was found that, 
although the inputs were not purely sinusoidal, the steady state condition 
was rapidly reached. Harmonic analysis of both the responses and the 
forcing function was considered but was discarded as being too time con- 
suming. Since all the data had been recorded on magnetic tape, a tech- 
nique of filtering both the forcing function and the aircraft responses 
with a low frequency band-pass filter was developed. The procedure, as 
previously described, was really nothing more than a fast and simple 
method of performing harmonic analysis of the traces and required that 
the data be recorded on magnetic tape. The technique was a valuable 
one, and it was felt by the authors that it could be utilized . to advantage 
in all forms of steady state dynamic flight testing where telemetering 
and magnetic recording of data is performed. 

The use of telemetering equipment in a flight test program of this 


sort requires more complicated instrumentation than the simple SFIM 
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recorder utilized in Ref. 1, in that a ground station and its associated 
equipment is required. The actual instrumentation in the aircraft is 
little complicated by the substitution of the telemeter transmitter for 
the in-flight recorder. The utility derived from having magnetic tape 
recordings when the telemeter scheme is used is enormous, and it is 
felt that the advantages of this method outweigh the complications pre- 
sented by the extra instrumentation required. The filter technique of 
course would be impossible without tape-recorded data, but in addition 
tape records of this sort may be expanded to any desired time scale 
making phase angle measurements easier and more accurate. It is 
also possible to amplify small responses to distinguishable values, but 
care must be taken not to exceed the accuracy of the record. The tele- 
metering scheme utilized here was accurate to 1% of full scale and thus 
small responses could be expected to have a much greater percentage 
error than large ones. It was therefore desirable during this investiga- 
tion to obtain responses which were as large as possible without exceed- 
ing the range of linearity of the equations of motion. This was accom- 
plished by varying the forcing function amplitudes over the frequency 
range to give desirable responses. In general, large control deflections 
were required at high frequencies and small ones at and below the dutch 
roll frequency. 

The method used in this investigation to determine stability deriva- 


tives from frequency response was basically sound. As mentioned 
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previously, considerable difficulty was encountered in the determination 
of the cross derivatives, and it was felt that accurate determination of 
these would require much more accurate frequency response data than 
could be obtained with the simple instrumentation utilized here. The 
accurate determination of phase angles from data obtained using the 
forced sinusoidal oscillation technique is difficult indeed, and phase angle 
errors can have a large effect in the equations for the solution of the 
stability derivatives, especially when the coefficients of these equations 


involve sines of angles near 0° or 180° or cosines of angles near 90° 


or 270°. The technique developed here, which involved fairing a 
smooth curve through the experimental points and taking data for the 
computation of the derivatives from the smooth curve, eliminated some 
of this difficulty, although it was felt that the method of averages could 
probably eliminate most of this random error even if individual data 
points were used. 

The equations derived from the rudder data were, as previously 
mentioned, ill-conditioned for solution, and this fact required determin- 
ation of some of the derivatives from the aileron data before the rudder 
equations could be used. Intheory, it is possible to determine all the 
derivatives, except of course the aileron control derivatives, from 
rudder data alone. While this was not possible for the Navion aircraft 
at the condition investigated, it was felt that for other aircraft or other 


flight conditions the difficulties encountered here might not be present. 
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It became apparent to the authors that a need exists for extensive study 
of the determination lateral stability derivatives from frequency re- 
sponse data, especially as concerns ill-conditioned equations and possi- 
ble methods for their solution. While it finally proved possible to 
obtain fairly accurate derivatives in this investigation, certain devia- 
tions from the basic method were required. The methods applied here 
were peculiar to this particular aircraft and flight condition and might 
very well prove unsuccessful for other aircraft or flight conditions, The 
best procedure to follow in the extraction of stability derivatives from 
frequency response data is, at the present time, a function of the type 
and accuracy of the data, although the authors feel that with time a 
standardized procedure for eliminating difficulties could be developed. 

The use of the method of averages for eliminating random error 
instead of the method of least squares, as outlined in Ref. 4, proved 
quite satisfactory and had the added advantage of reducing the computa- 
tion time considerably. When using this method care must be taken 
when grouping the points for averaging. If automatic computing facili- 
ties are available and large amounts of data require reduction, the 
method of least squares would probably prove more satisfactory and 
less subject to error. 

It was difficult in this analysis to acquire a feel for the accuracy 
of many of the individual derivatives. As previously stated, the nature 


of the equations and the fact that small differences of large quantities 
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were involved in the solution for C Ci,» and uS made the 


Ya ' 
accuracy of these derivatives questionable. In the solution for Cy; > 


Cir. , and C the controlling terms involved “ሃሪ C1, and 


n (a 
Cay respectively, and thus these control derivatives were also of 
doubtful accuracy. Lateral acceleration data would have eliminated 

this difficulty in the case of the side force derivatives, but was not 
obtained here. In the determination of the other derivatives, however, 
the potential degree of accuracy was not evident from the equations. 

A solution of the characteristic equation using the experimentally 
derived derivatives resulted in a period of 2.38 seconds and a time to 
half amplitude of 1.69 seconds for the dutch roll mode. Transient data 
obtained for the same aircraft at the same flight conditions indicated a 
period of 2.40 seconds and a time to half amplitude of 1.70 seconds. 

In view of this close agreement, it was concluded that Cn g , which 
primarily determines the period of the dutch roll mode, and Cnr’ the 
primary damping derivative, were accurately determined. 

A comparison of the frequency responses derived from the experi- 
mentally determined derivatives and the experimental points indicated 
certain regions of discrepancy. In Fig. 19 the calculated amplitude 
response in roll to aileron oscillations was noticeably low at and above 
the dutch roll frequency, and the calculated phase angle was slightly 
high in the same region. An analysis of the transfer function contained 


in Table V indicated that an increase in the coefficient of the first power 
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term in the numerator would remedy this condition. The most obvious 
means of accomplishing this was an increase in the value of Cis, 

An increase in this derivative would also increase the constant term in 
the numerator of the yaw transfer function and thus increase the ampli- 
tude response at low frequencies. Fig. 20 shows the need for this cor- 
rection. While a detailed analysis of the effect of increasing Cisa 

was not performed, these considerations led to the conclusion that the 
value determined for 61 was low, and a comparison with the results 
of Ref. 5 indicated this to be the case. 

In Fig. 22 the calculated amplitude response in roll to rudder 
oscillations is substantially less than the experimental response at low 
frequencies and the phase angle is consistently larger for all frequencies. 
An analysis of the transfer function for this response indicated that an 
increase in Cn, P would tend to correct both these discrepancies. In 
addition, the calculated low frequency amplitude response and high fre- 
quency phase response for sideslip, as shown in Fig. 24, would also be 
altered favorably by a larger Cn, . A comparison with the results 
of Ref. 5 indicated that the value determined here for Cn gy Was too 
low. 

Since the match between the calculated and actual responses was 
generally good with the exception of the cases mentioned above, it was 
felt that the determination of those derivatives with the greatest effect 


on the frequency response, i.e., رت‎ , <a, > ole and Cn, : was 
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good. While detailed curve fitting of individual responses by altering - 
individual derivatives was not considered within the scope of this inves- 
tigation, such a technique for correcting derivatives might warrant 


further study. 
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CONCLUSIONS AND RECOMMENDATIONS 


The method utilized in this investigation to determine the lateral 
frequency response of an airplane proved entirely feasible and could 
undoubtedly be applied to many current production aircraft. The 
method, however, had the disadvantage of requiring more flight time 
than currently popular transient response methods, and thus would not 
be practical for missile testing. 

Accurate determination of frequency response previously required 
an elaborate system for the generation of an exact sinusoidal forcing 
function or time consuming harmonic analysis of the data. The filtering 
technique developed here eliminated these difficulties, and could have 
wide application in this type of flight testing. 

The method utilized here for the determination of stability deriva- 
tives from frequency response data proved accurate for those derivatives 
with significant effect on the frequency response. It was, however, im- 
possible to determine Cn, * C1, : Cn, , and “14. with a high 
degree of accuracy. Better results could have been realized in the 
determination of “ሃሪ and Cy, if lateral acceleration data had been 
obtained, and it is recommended that such measurements be made in 
future tests. It is further recommended that a study be made to develop 
a systematic method for eliminating the difficulties encountered here 


with ill-conditioned equations. 
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TABIE I 


PHYSICAL CHARACTERISTICS OF NAVION 


Ae WING 


1. Total Wing Area, Sy (includes flaps, ailerons 
and 19.87 ጀቲት covered by the fuselage) 


2 Span By 
3. MAC 


Angle of Incidence‏ وبا 
Root, ir‏ 
Tip, it‏ 


5. Twist 
Aerodynamic 
Ge ome tric 

6. Airfoil Section 
Root 
Tip 

7۲۰ Aspect Ratio, AR, 

8. Taper Ratio, A y 

9. Dihedral 

10. Root Chord 


11. Tip Chord 


A-l. AIIERON (one Aileron only) 
(Beaded skin, Frise Nose Balance, 


No Trailing Edge Extension) 
1e Area, یت‎ 

2, Span, ba 

3. Aileron Deflection, G a 
Le Boost 


So Trim Tab 


181.337 Ft.2 
33.378 Ft. 
68.352 in. 


+20 
-10 


29 31! 
39 


NACA ۶ 
NACA 6410R 
6 Ol 
0.5265 

7.5° 
7.2 Ft. 
3.92 Ft. 


1200-1200 


Wheel Throw 
Fixed bend tab on 
right aileron. 


2,161 Ft.? 
61.987 in. 
309 UP, 209 DN. 
None 


Fixed bend tab on 
right aileron. 








Bo 


8-1 0 


B-2 . 


C. 


HORIZONTAL TAIL 


Total Area‏ ول 
(Includes 2.368 Ft.2 covered by fuselage)‏ 


2. Span, by 

3: Airfoil Section 
4. Taper Ratio, Ay 
5. MAC, Ky 


6. Aspect Ratio, R t 


HORIZONTAL STABILIZER 


1, Area, Ss 


2. S u (with reference to FRP) 


ELEVATORS 


(No Trailing Edge Extensions) 


lo. Total Area, Se 

2. Span, be 

3. Deflection, De 

Trim Tabs (32" Span, h-1/2" Chard)‏ وبا 
Root Chard‏ .5 


6. Tip Chord 


VERTICAL TAIL 


l. Total Area, Sy 


(Includes 2.577 Ft.? blanketed by 
fuselage and excluding 1.183 Ft.2 
af dorsal fin) 


2. Airfoil Section 
Root 
Tip 


13.051 Ft.* 


13.172 Ft. 
NACA 0012 
.67 
3.3 Ft. 
با‎ 2 


28.953 Ft.? 
0ه‎ 

Smooth Skin 

Flat sided 

No trim bungee 

Balance spring 

1.098 Ft.® 

73.582 in. 

300 UP, 209 DN. 

300 UP, 300 DN. 
1.5 Ft. 


1.0 Ft. 


12,925 Ft.? 


NACA 0013.2 Mod. 
NACA 0012-6), " 





0-1. 


6-2 ۰ 


5. 


Eo 


0. 


VERTICAL STABILIZER 


le Area, Sy, 6.873 Ft.2 
(Includes 0.1427 Ft.? blanketed by fuselage) 

2. X (With reference to fuselage f) 20 Nose left 

RUDDER Smooth Skin 
Rigged 3° Rt. to 
Fin f, Fixed 
Bend Tab. 

le Area, Sr 6.052 Ft.? 

2. Rudder Deflection, © r 17° L - 239 R 

3. Rudder Pedal Throw 5.75 in. 

Lo Trim Tab Fixed Bend Tab 

PROPELLER Hartzell 

lo Activity Factor, A.F. 100 

2. Diameter, d 86" 

3. Pitch 00.75R 21° 301 

NOSE GEAR 

l. Travel + 30° 

2. Strut Angle (To Vertical) 109 

MISCELLANEOUS 

le Length, overall 27.25 Ft. 

2. Tail length, AE 16.88 Ft. 

POWER 

l. Engine Continental E-185 Engine 


185 HP @ 2300 RPM, 29" Hg 
M.P. @ Sea level 





TABLE II 


EXPERIMENTAL FREQUENCY RESPONSE DATA 


TO SINUSOIDAL AIIBRON OSCILLATIONS 


۳۰۸ ۰5 م‎ e 129.6 m.p.h. 
Density Alt. = 5000 ft. 
Gross Wt. = 2870 lbs. 


o IRI | ኢ|፪) | *..| ISL | ee 
0.73 ۱ 3,20 9 .0 | 1.225 | - 11.0 | .362 173.3 
1.14 | 2.175 | 76.0 .703 | 7.5 | oh & | 152.2 
1.61 13,552 | 88.2 | 637 03 28.9 | 566 ¦ 110.0 
1.78 | 1.322 | 85 .0 | «60 30.0 | .529 | 145.0 
223 | id] | 3 | .ዜዜ | -ቧ.ቹ0 | .ጩዐ | 108.1 
2.33 | 1.012 | 8.0 „881 ۱ - 58.0 | „806 | 119,0 
Brem | 1.197 | alo 790 | -308.0 | 1.032 | 72.8 
2059 | 1.458 | 88.5 | 750 | -132.2 | 1.000 46.0 
2.62 | 1.130 | باء ب101‎ | „820 | - 91.0 | 1.100 | 85.0 
2.79 | 148 | 900 | .T5 | -136.7 ,ፆቿ | 16.0 
2.97 | 1.338 | 83.2 | ,66.  -332.8 | 85 | 55. 
3.27 | 1.001 | 69 باء‎ | 2259 | «180.0 | — .367 | 8.7 
319  .89 | 66.1 | .23 | 172,0 283 | 10.0 
3.92 | .822 | 65.3  .1608 | — 166.6 .260 | 356.5 
1.63 | .681 63.5 | 40853 151.0 — E | 2 
5.20 |  .562 77 13.3 | .1265 | 33.7 
5.85 | 96 ' 89.0 | .0l08 125.0 | .0775 | 336.3 





6 were 9 
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فو > سے سر ای مہ یھت ہے سو با eee Cee en eee ete Semen‏ سرب و uon‏ ہے .2.0 - 


TABLE III 
EXPERIMENTAL FREQUENCY RESPONSE DATA 
TO SINUSOIDAL RUDDER OSCILLATIONS 
Te... = 129.6 10۰0912۰ 


Density Alt. = 5000 Ir. 
Gross Wt. = 2870 lbs. 


| a e, | da | 5766 
E ور‎ Br اھ‎ $e ۱ ا‎ M Leg 
1.668 | 1.222 1002 7.1 | 176.2 | 1.175 6.1 
| | | | | 
31:235 ጃ 14121 94.5 . - 90h | ፲68.ከ 1,130 | - ከ. 
| 1.615 : 1.100 94.6 1.053 | 171.8 irs ۰ 


1.8 ` 1.4106 7904 1.211 : 160.0 
82.5 1.341 | 
2.21 1.358 $71 ! 1.810 | lo | 2.135 ا‎ 
2.1 1065 37.3 | 2.165 ! 339.9 ۸۱ 2.55 
248 — 13,8 . ጪሪ | 2.3. | 132.5 | 2.60 


2.58 1.55 | 21.6 | 2.21 | 106 ¡ 2.63 | 29908 


TON nn 


1.87 ; 3,352 
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2.5. . 3.3. | 35 ; 2168 | 124.0 Í 26 | 308.6 
2.69 1.547 +5,2 0 | 79.5 | 5 | 281.1, 
2.99 1.00 | _2he5 ^ 10.661 . 2 | 1.99 | 216.2 
aed 452 | -51.8 | 822 1. 38.9 | 980 | 221.9 


۱:00 .29b | -63.6 | .599 | 29.2 ¡¦ ,ይ55 | 210, 3 
وود‎ ' 32 | -፳? | 48 ! 20.0 | ,3፳ | 
11.68 .181 -66.6 0390 | 30.9 ` ,390 | 226.5 
TS ag  ,[ዐፀ | 215.7 
lı. `  .160 | 68.0 .338 | 200% - E 


5.30 38 ፤ -776 . ,5 . 36 ! سرت‎ በ5 


$38 5 =- ! - ! 235 22.6 | .85 i ሟጨ8 
5.73 | .083 | -7505 | .216 | 3ከ,5 | . 3 ፡ 208,7 
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i à : 
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TABLE IV 


INITIAL THEORETICAL AND EXPERIMENTALLY DERIVED 


LATERAL STABILITY DERIVATIVES FCR THE NAVION 


TE Ador E 129.6 MePehe 
Density Alt. = 5000 ft. 
Gross Wt. « 2870 lbs. 


Initial 
Theoretical 


Derivative 


-0.592 


+0. 102 


“0, 083. 


-0.160 


-0.0937 


+005 


-0.051 


+0.123 


-0.101 


+0.0152 


+0.0138 


-0.0625 





Results 
of Ref. 5 


-0.592 


+0. 131 


-0.0595 


-0 0 


-0.0951 


-0.0705 


بابلا .0+ 


-0.137 


+0.0069 


+0:0031 


0.077 





-0.0513 


-0.13l, 


-0.0908 


+0.0895 


-0.106 


+06 


-0.108 


+0.0128 


+0 .00376 


-0.0673 


RS SR ae ee A AA‏ ی 





(3) 


(5) 


(3) 


(ከ) 


(5) 


(6) 


ری فک 
Density Alt. = 5000 ft.‏ 


Gross Wt. 


2 
da 


مج 


XE 


51m 


وه 
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TABLE V 


TRANSFER FUNCTIONS FOR THE NAVION DETERMINED 


FROM EXPER IMENTAL FREQUENCY RESPONSE DATA 


= 129.6 m.poh. 


= 2870 lbs. 


-.244 5” - 20.0415 - 2.155 
s + 


کہ ہے نو مس መዌጨ ፍው‏ ی 


+ 9.002 5° + 13.568 57% + 57.8645 -- ۱ 


- 23,470 5* — {7.2355 — 136.599 


55 + ዓ.ዕ6255 4 |3.ናዊሯ 5” ¿5764 35-103 
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S (54+ 3.00259 «13.568 s^ - S1, $64 S - 4631) 


S^ ı 41.9915 -1.450‏ 11.96% + ہی 
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SUS 4 9.050239 »/3.5649 ^ ^ 7. Pf9MS—^/^. 63l 
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APPENDIX A 


Development of Vector Equations of Lateral Motion 


for Determination of Stability Derivatives 


The three linear equations of Lateral Motion may be written in 


the following form: 
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Roll: (2) Ga gu ርፉ ) D — = ۶ D*+ ርረ D) > CH ác ረር ፡ $a‏ 
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The coefficients of these equations are defined as follows: 
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The vector technique is applied to the equations in the following 
manner. If an input Se to a linear system is considered to have a 
sinusoidal variation of frequency ù) , it can be considered to be a 
vector of unit magnitude lying along the real axis of the complex plane 


as shown in the sketch below. 
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The first derivative of the sinusoidal input is obtained by multi- 
plying amplitude by W and rotating the vector through 90°. The 
steady state response of the system may be considered as a vector of 
amplitude ratio R and having a direction displaced from the real axis 
by a phase angle, $ . Thus, in the complex plane, the vector for 


sideslip may represented as 


Dc |‏ ہے کا 
B= ka e * 4 E (cos $e, + 4 Sin $2, )= Agt+ Be‏ 





In complex notation the derivative of either the control input or the 


response for the steady state condition may be determined by multiply- 


ing by «wu and the integral by dividing by Au à 


The following are defined: 
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In this notation with the input vector having a magnitude of unity the 


responses and their derivatives may be defined as follows: 
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Substitution of equations (4) to (12) in equations (1), (2), and (3) and 
separation into real and imaginary parts yields the following six 


simplified equations: 
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In these equations the A's and B s will be available from the 
frequency response curves for any chosen value of w . In theory it 
is seen that equation (13) may be solved for er from the response 
at a single frequency since Cue is known. Equation (14) may then 
be solved for Cy. using the computed value Cyz * The complete 
solution of the yaw and roll equations may be accomplished if the 
responses are known at three different frequencies since the maximum 
number of unknown coefficients in any one equation never exceeds 
three. The method of averages attempts to eliminate random error in 
single points by writing each equation for a number of frequencies and 
grouping the equations thus determined in as many groups as there are 
unknowns. The coefficients for each group are summed and the result- 


ing simultaneous equations are solved for the unknowns. 


In the roll and yaw equations it is found that better solutions are 





obtained ifthe grouping of points is performed in such a manner that 
all the points in a given group lie in the same frequency range and 
further that the three groups chosen include data below, at, and above 
the dutch roll frequency. 


The final computational equations then become: 
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The development of the computational equations was performed 
۳ * 
for the rudder case with the coefficients Cy رم ی‎ and Cn 5 equal to 
D ረኒ- 
zero. It is obvious that for the aileron case the equations are identical 


ታ 
to those developed above, with the substitution of Cy, and Cp 7 
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for Cif and NÉS . 
or C], an Cas. ጋ ከ ከክ... 


If difficulty in the solution of the equations arises because of ill- 
conditioning, it may be possible to resolve in some direction other than 
that of the control deflection vector. It may also be possible to elimin- 
ate some troublesome coefficient between equations and thus obtain 
solutions. No attempt was made here to consider all the possibilities 
of modification of the equations. 

If the derivatives are known, the transfer function coefficients 


may be easily developed as follows: 
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follows, where 


The coefficients of the transfer functions are determined as 
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Once the transfer function coefficients are known the frequency 


response may be determined by the substitution of Cw for s, giving 


the following: 
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